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Einstein-dilaton-Gauss-Bonnet gravity is investigated on existence of solutions with mild sin- 
gularities, not shielded by the event horizons. These still may have sense since presumably such 
singularities will be smoothed by corrections to Einstein theory from quantum gravity /string theory. 
We show that gravity with the first-order correction, the Gauss-Bonnet term, gives rise to special 
types of singularities, which we call 'cusps', with 1/2-th and 1/3-th powers in series expansion. The 
full space then can be split onto several cusps with classically impenetrable borders, and/or flat 
asymptotic. 



^ , T. GAUSS BONNET GRAVITY 

^ ■ Gauss-Bonnet term provides a non-trivial contribution to the Einstein equations only in dimensions D > 4, while 
in D = 4 gravity one should include interaction with the dilaton field. The solution without the GB correction is 
a standard Schwarzschild metrics with constant dilaton (zero dilaton charge is ensured by the 'no- hair' theorem). 
Adding the GB term gives rise to a non-trivial dilaton configuration with two different types of asymptotically flat 
solutions: the first is the black hole, and the second contains a naked singularity. Actually the BH solution also 
_^ ■ contains singularity, but this is shielded by the event horizon. 

Solutions with naked singularities are usually not in a favor, therefore subsequent investigations were produced only 
for the black hole branch. But now let's turn back to our motivation of adding a higher order curvature corrections to 
Oh' the action. It is because of our belief in some fundamental quantum theory we add the GB term. Therefore when we 
get a 'singularity' this just means that we moved outside the semi-classical limit with first-order corrections, but the 
'real' quantum theory should deal with the growing curvature and energy density in some way. This is why the domain 
of initial conditions with high energy density should not be forgotten in favor of the domain of initial conditions with 
the event horizon. So in our work we will explore the family of solutions with GB corrections in more details. 

Again, one should mention that since the GB term is motivated by string theory, the corresponding calculations 
should be produced rather in the string frame than int the Einstein frame, and the transformation of the Gauss- 
Bonnet term between these frames generates a very complicated additional term. It was shown that not all solutions 
■ known in the Einstein frame are reproduced in the string frame the string frame [l| in the context of the dilaton black 
holes with Gauss-Bonnet corrections. For example, the family of singular solutions, obtained by (H as far as we know 
was not found in the string frame. 

Therefore in our investigation we will provide two parallel calculations: for the usual Einstein-dilaton-Gauss-Bonnet 
action (to extend known results) 
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S [E) = ^ J (R - (0 M In S)7(2a 2 ) + aSK 2 GB ) y/=g d 4 x, (1) 

and for the stringy version of the action 

S {str) = ^-J(R+ (d. In Sf/a 2 + aH GB ) S^9d A x. (2) 

Note that in [l[ the first version was called a 'truncated' one. We shall use another notation: the first version will 
be marked as "(E)" (EGBD), and the second version as "(str)" (SEGBD). Along with the Ricci scalar curvature R, 
it contains the dilaton field S = e 2a ^ and the Gauss-Bonnet term 1Z GB . The Euler density of the GB term has the 
following form: 

T^GB — R 2 ~ ^RpvR 1 * 1 ' + Ra/3fivR al3fil/ ■ (3) 

The action is written in Plank units with h = c = G = 1, and it contains two parameters: a which is a dilaton coupling 
constant and a which is a 'correction' parameter in sense that GB term is treated as a correction to the Einstein's 
action. Both parameters are supposed to be positive. 
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The static spherically symmetric spacetime can be described by the following metrics: 

dr 2 

ds 2 = -w(r)a(r) 2 dt 2 + — - + p(r) 2 dn 2 (4) 
w(r) 

which contains three functions depending only on the radial coordinate, but one of them can be gauged away by the 
rescaling of r. 

One can transform the action in the string frame ^ to the Einstein's frame by the conformal transformation 



S ( q str) = / R\g$S- 1 ]Sy/-gV s) S- 4 d*x 



- 3/2 / In S) 2 V Z 9^ d A x. (5) 



In string theory one has a = 1, therefore Sg Str ^ + sjf — Sg + si , where Sd is a dilaton kinetic term. So in the 
case a = 1 the standard Einstein-dilaton action without GB corrections is identical in both frames up to the conformal 
transformation. For arbitrary a the transformation from the string frame to the Einstein frame produces the action 
with additional kinetic term. 

The transformation of the GB term under the conformal transformations is rather complicated. But we can just 
write the GB-term in the non-transformed frame and substitute the metric functions as w 

(str) = Sw (E) 

S^V^, p( 8tr ) = S-V 2 p( E \ This will allow us to express both versions of the action in one frame, where the 
high-order curvature correction terms will be different. 
Since the GB term is a total derivative: 

^Tl GB = A' GB , where A GB = ~(wa 2 )'(wp' 2 - 1), (6) 

(7 

it enters the lagrangian like —aS'A GB . The expression above is correct for the gju> metric functions in the string 
frame and for the gjfjy metric functions in the Einstein frame. Therefore the GB part of the action l S < - str - ) in the 
Einstein frame will enter also like —aS'A GB , but with the different expression for the Acs-term which now reads as 

(7) 
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where the metric functions were substituted by the transformed ones. This expression can be split onto the same 
Acs-term when the derivatives do not act on the dilatonic function and the terms which are proportional to it's 
derivatives: 

4 

A GB -A GB + AA GS , where AA GB - £ A„(ln S)' n . (8) 

n=2 

One can see that the difference between two forms of the GB term is significant only when dilaton demonstrates the 
exponential growth, and it vanishes if asymptotically S — >■ const. 
Now we can obtain the one-dimensional effective EDGB lagrangian: 



CT (- 2 )W 2 -1). (9) 



r (E) p'{o 2 wp)' awp 2 S' 2 aS' 2 , 

C ~— 8^"— (W) 

And the lagrangian of the stringy model reads as 

r( *r) = s ( p^pL + « + ^s^\ + ^vr _ ^ (w2)W2 _ 1} . (10) 

\ 2(7 2 4cri> z / 4ap z a 

The first solution for the Einstein DGB system was obtained by Kanti et. al. Q where they obtained the BH-branch 
of solutions with horizon and another branch with an unshielded singularity r x where square of the Riemann tensor 
diverges like (r— r^) -1 . Then, Alexeyev et. al. [H explored the behavior of the solutions under horizon Th and revealed 
that then r s point becomes a "turning point" singularity from which two branches of solutions arise. The first one 
goes outward from r s to r/j and the second one moves further into the black hole from r s to the singularity r x . In 
[J] the dynamical analysis was produced and it was shown that the only regular solutions with flat asymptotic are 
the black holes. In paper [l[ the authors reproduced the BH solutions in the string frame and revealed the difference 
between EDGB and SEDGB models, but they did not investigate non-BH solutions in details, just mention that 
sometimes there is not a BH, but a singularity. 
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II. CUSP SINGULARITIES 



We treat cusp as a mild singularity of the metrics, where the dilaton function and metric components are regular 
and non- vanishing, but their second derivatives are singular. 

In what follows we will use the shifted radial variable x = r — r s for convenience. In the vicinity of the singularity 
the metric and dilaton functions can be written as expansions by powers of x: 

oo oo oo 

™ = E W «A X " /Z ' P = J2pn/^ n/Z , S = ]>>„/ 2 *" /Z , (11) 

n— n=0 n— 

with z = 2, 3, . . . parameterizing the order of singularity. First, we will investigate the case z = 2, which corresponds 
to the turning point of the metrics. To provide the regularity of the first derivatives we impose the condition 
w i/2 = P1/2 = S1/2 = 0. In this case the curvature has 1/v^ singularity at the cusp. 

Substituting the expansions in to the equations of motion we find that starting from the second order (i.e. n/2 > 2) 
the equations are linear by the n/2-th coefficients, expressing them through the previous coefficients. Unfortunately, 
in the first order (to be precious 3/2 order) we obtained the system for the P3/2 and 93/2 which contains polynomials 
of orders higher than five. Therefore we are unable to provide explicit expansions. Nevertheless, this system can 
be solved numerically with any required precision when substituting quantities instead of free parameters of the 
expansions. 

The free parameters of the expansions are the parameters of the configuration (a, a), the dilaton value on the cusp 
surface, pa, the radius of the cusp, pa, the angular deficit of the cusp which can be expressed by p\ and the time 
component of the metrics, wq. But the last two parameters can be combined into one u = wap\, while the remaining 
degree of freedom will be responsible only for the scale transformations. In the case of flat asymptotical metrics the 
GB term vanishes faster then the curvature term, for which it is well known that the asymptotic should satisfy the 
relation Woop'^ — 1, while separately and p'^ can be not equal to unity. Therefore the parameter u can be treated 
as a deviance of the cusp metrics from the flat asymptotic. Also the dilaton parameter p can be excluded by the 
appropriate rescaling of the dilaton function and the GB parameter a — » apa- In what follow we will set pa — 1, so 
that the effective GB parameter apo will coincide with the initial parameter a in the lagrangian. 



A. z=2 cusp 



The most interesting solutions are those which interpolate between the cusp and the flat asymptotic (Fig. [T]), and 
numerical computations confirmed their existence. Other solutions interpolate between two cusps of different radii 
(Fig. [2]). The actual parameter space is four-dimensional (a, a, pa, u) but it can be schematically described in the 
following way. The key role in producing the asymptotically flat solutions belongs to the parameter u. It appeared 
that the deviance from the asymptotic should be small enough, but non-vanishing. There is a forbidden gap around 
u = 1, and two permitted bands on each side of it. The decrease of a, a and the increase of pa causes the stretching of 
the allowed bands, while the forbidden gap also enlarges above u = 1 and slightly modifies below u = 1. Conversely, 
the increase of a, a and the decrease of po causes the contraction of the both permitted bands and the forbidden gap 
above u = 1. To be precise, these bands are not solid, they possess stripy structure with narrow forbidden gaps. But 
the detailing seems to be abundant for practical purposes. 

One can conclude that all conditions which stimulate existence of asymptotically flat solutions can be understood 
as relative decrease of the contribution of Gauss-Bonnet corrections in comparison with the classical Einstein-dilaton 
gravity terms. This can be seen from the table which follow. Data in each table starts from the point a = a = pa = 1 
and then one of the parameters varies in the direction which enlarges the allowed bands for u. 



(a, a, pa) 


allowed band 


allowed band 


(1, 1, 1) 
(0.5, 1, 1) 
(0.33, 1, 1) 
(1, 0.5, 1) 
(1, 0.33, 1) 
(1,1,2) 
(1, 1, 3) 


0.93 0.99 
0.77 0.99 
0.65 0.98 
0.89 0.98 
0.75 0.94 
0.46 0.92 
0.6 0.86 


1.02 1.11 
1.09 1.22 
1.1 1.34 
1.06 1.23 
1.11 1.34 
1.16 1.55 
1.54 2.5 



One can compute the ADM mass and the dilaton charge of the configuration using the asymptotic obtained nu- 
merically. It appears that the dilaton charge decreases and the ADM mass increases when the parameter u deviates 



4 



from unity. The dilaton charge is negative, while the mass is positive on the lower band and negative for on upper 
band. 

There is another interesting branch of solutions for the considerable deviance from the asymptotic with u ~ 0.6. 
With growing x one obtains solutions with the metrics reaching the flat asymptotic, while the dilaton function grows 
exponentially. Still these solutions are singular, but the point of singularity moves to infinity with the decrease of 
a. Actually the plot of the dilaton function looks like that of the scale factor in the inflationary scenario with fast 
acceleration period turning then into slow deceleration. With the significant decrease of a this branch disappears 
since for sufficiently small a the lower allowed band of regular solutions reaches the point u ~ 0.6. 

B. z=3 cusp 

One can choose the value z — 3, i.e. the expansions near cusp will contain powers of n/3. In this case the condition 
of regularity of the first derivatives does not hold, but there is a w 2 /3X 1 ^ 3 term. The scalar curvature possess stronger 
singularity a; -4 / 3 instead of the square root singularity as it was for z = 2. 

Unlike the previous case, the coefficients of the expansions can be found explicitly, but they are rather complicated. 
It turned out that there are two branches of the expansions depending on whether P4/3 vanishes or not. Moreover, 
each branch splits into two other branches which are identical up to the transformation x n l z — > ( — l)"x"/ 2 , which 
defines the 'direction' where the solution spreads. The parameter space contains five parameters, but again po and 
wq can be scaled away. So the effective parameter space (a, a, po) is three-dimensional, and there is no independent 
parameter like u in previous case. 

Surprisingly, we did not found any regular solution with flat asymptotic numerically. Each solution ends by the 
square-root cusp considered in previous section as shown on (Fig. [3]). Nevertheless, there exist the same branch of 
solutions with the exponentially growing dilaton function when a decreases. But now the absence of regular solutions 
does not bound this branch from below. So with a —¥ po on has S — > oo and w, p' — > const what gives us the flat 
metrics. The point of singularity moves to infinity. Mention that there is no cusp for a = 0, since the terms with 
non-integer powers of x vanish in the expansions both for z = 2 and z = 3. 

III. CONCLUSION 

We investigated the class of singular solutions to Gauss-Bonnet dilaton gravity system. First, we mentioned already 
known singular solution with I/2-th powers expansion near singularity. Such solutions can provide flat asymptotic, 
but also they can demonstrate from-cusp-to-cusp behavior, where the entire space can be split onto several regions by 
cusps which are classically impenetrable. Then we showed the existence of a new singular solution to Gauss-Bonnet 
dilatonic gravity, containing 1/3-th powers near the singular point. It can not produce flat asymptotic and always 
spreads from one singularity to another. But it improves the inner region in the sense that one has flat space between 
cusps. 

Still there is an open question if there are other cusps, probably even with not a power-series expansions. Such 
nonlinear differential equations can provide a very diverse picture. We investigated 1/4- and 1/5-th powers and found 
no cusp solutions. Another question which we explored was the difference between classical and string versions of 
dilatonic gravity with curvature corrections. It appeared that qualitatively the picture is almost the same. 

One can ask about applications of cusp solutions. Some of them can be interesting in view of the brane-world 
scenarios, where the singularity is moved to the bulk 0-0] ■ Finally, if the black holes were treated as some unnatural 
solutions for many years, and now we have several dozens of candidates to them, probably another even more strange 
configurations will arise in future astrophysics. 
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